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We present near-infrared imaging obtained with ESO VLT/ISAAC of a sam- 
ple of 16 low luminosity radio-quiet quasars at the epoch around the peak of the 
quasar activity (2 < 2; < 3), aimed at investigating their host galaxies. For 11 
quasars, we are able to detect the host galaxies and derive their properties, while 
for the other five quasars, upper limits to the host luminosity are estimated. The 
luminosities of the host galaxies of radio-quiet quasars at high redshift are in the 
range of those of massive inactive elliptical galaxies. 

This work complements our previous systematic study of quasar hosts aimed 
to trace the cosmological luminosity evolution of the host galaxies up to z ~ 2 
and extends our pilot study of a few luminous quasars at z > 2. The luminosity 
trend with cosmic epoch resembles that observed for massive inactive galaxies, 
suggesting a similar star formation history. In particular, both quasar host galax- 
ies and massive inactive galaxies appear mostly assembled already at the peak 
age of the quasar activity. This result is of key importance for testing the models 
of joint formation and evolution of galaxies and their active nuclei. 

Subject headings: Galaxies: active — Infrared: galaxies — Quasars:general — galax- 
ies: evolution 



1. Introduction 



There is nowadays compelling evidence that the processes of the formation and evolu- 
tion of galaxies and the nuclear activity are intimately linked. The most direct link at low 
redshift is the correlation of the mass of the central black holes with the luminosity (mass) 
and the stellar velocity dispersion of the spheroids where they reside, in both inactive and 
active galaxies (see Ferrarese 2006 for a review). If this link holds also at higher redshift, 
the observed population of high redshift quasars traces the existence of ^ 10^ M,^ super- 



massive BHs and m assive spheroids at very early (< 1 Gyr) cosmic epochs (IFan et al.l 12003 



Willott et al.l 120031). Moreover, the strong co smological evolution of the quasar population 
( iDunlop fc Peacocklll990l : IWarren et al.lll994j ) is similar to the star formation history in the 
Universe (Madau et al. 1998; see also Lapi et al., 2006). According to the hierarchical 



merging scenarios for structure formation and evolution (e.g.. iKauffmann fc Haehneltl 12000 



Di Matteo et al.ll2005l ). the massive spheroidals should be the products of successive merger 
events, and the properties of quasar hosts (mass, luminosity, size) should show a redshift 
dependence. Moreover, the quasar activity is expected to deposit large amounts of energy 
in the spheroid, possibly suppressing star formation (e.g.. Silk & Rees 1998). 



In this context, the direct detection and characterization of high redshift quasar host 



- 3 - 



galaxies is crucial to understand the joint cosmic assembly and evolution of bulges of galaxies 
and their central black holes. In particular, a key point is to probe the quasar host properties 
up to and beyond the epoch of the peak of quasar activity (2 < 2; < 3). 

Low redshift {z < 0.5) quasars are hosted in galaxies containing a luminous, massive 
bul ge component that becomes dominant in r a dio- loud objects and at high nuclear l uminos- 



itv (IHamilton et al.ll2002l : IPunlop et al.ll2003l : IPaeani et al.ll2003l : iFlovd et all 120041 1. Their 



stellar populations are believed to b e relatively old, especially in very lum inous AGN com 



Nolan et al. 


2001; 


Dunlop et al. 


2003) 



is some imaging and spectroscopic evidence for relatively young/intermediate age stellar 



(Kauffmann et al. 


2003 


: Jahnke et a 


. 2004 


Kotilainen & Falomo 


2004; 


Raimann et al. 


2005 


Hvvonen et al. 


2007bl: ] 


jetawe et al. 


2007 




Baldi & Capetti 


2008|; 


Hvvonen et al. 


2009). 



The detection and characterization of the host galaxies of high redshift quasars is chal- 
lenging since the quasar luminosity overwhelms the extended emission from the galaxy, 
especially in optical imaging, corresponding to rest-frame UV emission. Furthermore, the 
host galaxy surface brightness decreases rapidly with redshift. In order to cope with these 
difficulties, imaging with high spatial resolution and S/N together with a well defined point 
spread function (PSF) are essential. Systematic reliable studies of the host galaxies of z > 1 
quasars have thus became available only recently thanks to the use of near-infrared (NIR) 
imaging, where the nucleus-host luminosity ratio is more favourable, allowing to detect the 
old stellar population at high redshift. Concerning the stellar population, only few studies 
have been performed at high redshift (z > 0.5) so far. 



I n the largest study of quasar host galaxies available up to now (IFalomo et al.l 12001 



2OO4J : iKotilainen et al.l 120071 ). we carried out systematic imaging with the 8m Very Large 
Telescope (VLT) and the Infrared Spectrometer And Array Camera (hereafter, ISAAC) of 
32 quasars in the redshift range 1 < z < 2, to characterize their host galaxies. The sample 
includes both radio- loud quasars (RLQs) and radio-quiet quasars (RQQs), and it covers a 
wide range of the quasar luminosities. We found that the luminosity evolution of both RLQ 
and RQQ hosts is consistent with that of massive inactive ellipticals undergoing passive 
evolution. 

Beyond 2; ~ 2, detecting the host galaxies becomes extremely difficult, even with 



reported to 


DC resolved (Lehnert et a . 


1992; Hutchines et al. 19991: 


Ridewav et al. 2001; 


Croom et al. 


2004; Falomo et al. 2005; I 


^ene et al.ll2006; Falomo et al. 


2008; Schramm et al. 


20081: Villforth et al. 2008h mainlv folk 


jwing three different approaches: (1) observations 
rovide an excellent narrow PSF but are usually lim- 


from space fe.e.. iKukula et al. 2001). pi 
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ited by the small collecting area of the Hubble Space Telescope (HST). Furthermore, the 
PSF of the images taken with the HST WFPC2 are undersampled, leading to a system- 
atic overestimate of the flux from the host galaxy (see Kim et al. 2008). Note that this 
undersampling does not apply to images taken with ACS and NICMOS. (2) the extended 
emission from th e host galaxies is naturally magnified in gravitationally lensed quasars (e.g., 
Peng et al.ll2006l ). The drawback in this approach is that the PSF of the lensed targets is 
difficult to characterize, and the lens galaxy may contaminate the emission from the quasar 
host, making its detection extre r aely uncertain. (3) grou nd-based imaging with Adaptive 
Optics (AO) JCroom et all 120041 : IPalomo et allboosl . boosh usually satisfies the severe con- 



strains in the spatial resolution required to disentangle the extended emiss i on of the host 
galaxies from the nuclear one. In our previous studies (IFalomo et al.l l2005l . |2008| ). we ob- 
tained Ks-band images of quasars in the redshift range 2 < z < 3 using the AO system 
NACO at ESO VLT, which allowed us to clearly resolve two RLQs and two RQQs. The 
drawback in this approach is of statistical nature. In fact all current AO systems require a 
very bright star to be present close to the target. This therefore strongly limits the number 
of observable objects. This limitation should be overcome with the next generation of AO 
systems using artificial guide stars. (4) Another possibility, that we shall address in this 
paper, is to focus on high-z quasars of low luminosity where the nuclear-to-host luminos- 
ity ratio is more favorable. This ensures that, even without AO, the host galaxies can be 
detected and characterized from deep, NIR ground-based observations obtained in excellent 
natural seeing conditions. Moreover, the much larger field of view of non-AO images allows 
a significantly better characterization of the PSF from field stars. 

In this paper, we present a deep, high resolution imaging study with VLT/ISAAC of 
a sample of 2 < 2; < 3 radio-quiet quasars with relatively low nuclear luminosity, aimed at 
the study of their host galaxies. The structure of the paper is the following: In section 2 
we describe our sample. In section 3 we report the observations and data reduction and 
in section 4 we describe the data analysis. In section 5 the resulting luminosities of the 
quasar host are presented, together with a comparison to literature data available in the 
same redshift range. A discussion on the cosmic luminosity evolution of RQQ host galaxies 
is presented in section 6. Summary and directions for future work are given in section 7. We 
adopt the concordance cosmology with Hq = 70 km s~^ Mpc~^, Qm = 0.3 and = 0.7. 



2. The sample 



The sample of 2 < z < 3 quasars was extracted from the AGN catalogue of Veron-Cetty 
& Veron (2006) requiring: (a) rest-frame luminosity in the range: —26 > My > —27, k- 
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corrected assuming Francis et al. (1991) quasar SED templatq^, in order to maximize the 
likelihood for the host galaxy detection; (b) at least two bright (R < 16) stars within 40 
arcsec of the quasar, and many other fainter stars (R < 20, comparable to the brightness 
of the quasars) within the ISAAC 6.25 arcmin^ field of view, in order to provide a reliable 
characterization of the PSF. With these constrains, the search produced ~ 80 quasars, 
from which we extracted a sample of 16 RQQs with the most favourable observability from 
Paranal (Chile), covering the full redshift range. There is no statistically significant difference 
between the properties of the full and of the observed samples. Fig. [T] shows the distribution 
of the observed quasars in the z — My plane, compared with the samples of other similar 
studies available from literature and with the envelope of all the quasars in the AGN catalogue 
of Veron-Cetty & Veron (2006). 



Observations 



Deep Ks-band images of the quasars were obtained using ISAAC ( iMoorwood et al. 



19981 ). mounted on UTl (Antu) of VLT at the European Southern Observatory (ESO) in 
Paranal, Chile. ISAAC is equipped with a 1024 x 1024 pxl Hawaii Rockwell array, with a 
pixel scale of 0'.'148 px~^, giving a field of view of ~ 150 x 150 arcsec^. The observations were 
performed in service mode under photometric conditions in the period between 2006 March 
and September. The journal of observations is given in Table [H The seeing, as derived 
from the FWHM size of stars in each frame, was mostly excellent during the observations, 
ranging from ~0'.'4 to ~0'.'6 (average and median FWHM ~0'.'5, see Table [1]). Note that at 
the redshift of the targets, observations in the Ks-band probe the host galaxy at rest-frame 
5500 - 7000 A, roughly the R-band, thus allowing an easy comparison with optical studies 
at lower redshift. 

The in iages were secured using individual exposures of 2 minutes per frame, and a jitter 
procedure (ICuby et al.ll2000l ). which produces a set of frames at randomly offset telescope 
positions within a box of 10 x 10 arcsec. The total integration time was 38 minutes per 
target per observing block (OB). All targets but one (2QZJ124029-0010) were observed in 
two or three OBs of equal length. Since no significant discrepancy were found in the sky 
surface brightness and the seeing from the PSF of field stars between the different OBs for 
each target, we opted to combine all the individual images of each target. 

Data reduction was performed using our own improved version of the ESO pipeline for 
jitter imaging data (IDevillardl l200ll ) . Each frame was dark-subtracted and flat-fielded by a 



^Note that at 2 < z < 3 /c-correction is < 0.2 mag for a typical quasar SED. 
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normalized flat field obtained from twilight sky images. Sky subtraction was performed using 

median averaged and scaled sky frames obtained combining jittered exposures of the same 
field. Sky-subtracted images were aligned to sub-pixel accuracy, and co-added. Combined 
images arc trimmed into a 850 x 850 pixel (2.1 x 2.1 arcmin^) frame, covered by all the 
individual exposures. Finally, a polynomial surface was fitted to the combined image, after 
masking main sources, in order to remove spurious gradients in the background counts due 
to sky variations during the integration. 

Photometric calibration was performed from comparison with 2MASS magnitudes of 
bright stars available in the field. The estimated internal photometric accuracy is ~ 0.1 
mag. 

4. Two-dimensional data analysis 

Data analysis was carried out using AIDA (Astronomical Image Decomposition and 
Analysis; Uslenghi & Falomo, 2008), an IDL-based software package designed to perform 
two-dimensional model fitting of quasar images, providing simultaneous decomposition into 
nuclear and host galaxy components. The applied procedure is described in detail in Koti- 
lainen et al. (2007), and briefly summarized here. 

4.1. PSF modeling 

The most critical part of the analysis is the determination of the PSF model and the 
estimate of the background level around the target, which may strongly affect the faint signal 
from the object. To model the PSF shape, we used field stars in each frame, selected on 
the basis of their FWHM, sharpness, roundness and S/N ratio, including bright, slightly 
saturated stars, in order to properly model the faint wing of the PSF. A careful check of the 
light profiles and contour plot of each star let us exclude marginally resolved galaxies and 
stars with close companions. 

Each star was then modeled with four two-dimensional Gaussians, representing the core 
of the PSF, and an exponential feature, representing the extended wing of the PSF. Regions 
contaminated by nearby sources, saturated pixels and other defects affecting the images were 
masked out. The local background was computed in a circular annulus centered on the star, 
and its uncertainty was estimated from the standard deviation of the values computed in 
sectors of concentric sub-annuli included in this area. Finally, the region used in the fit 
was selected by defining an internal and an external radius of a circular area, allowing the 



- 7- 



exclusion of the core of bright, saturated stars. /^From the comparison of the resulting hght 
profiles (see Fig. [2]), no systematic dependence of the PSF was observed in the field of view. 
Thus, the same model was fitted simultaneously to all the usable stars of the image. 

The uncertainty of the PSF model was estimated by comparing the analytical fit with 
the individual observed star profiles, and adding a fixed term (0.1 mag/arcsec^) to account 
for possible systematic effects due to underlying assumptions in the data reduction (e.g. Zero 
Point and PSF stability, and perfect alignments in the de-jitter procedure). Panel e) in Fig. 
[2] shows the residual of PSF modeling together with the estimated uncertainty. 



4.2. Quasar host characterization 

In order to evaluate whether the targets are resolved, we first fitted each quasar image 
using only the PSF model. If the residuals revealed a significant excess over the PSF shape 
in the 0.4-1.0 arcsec range (where the contribution of the nuclear PSF quickly diminishes, 
while the contribution from the host galaxy signal becomes sufficiently high and detectable), 
the target was considered resolved. An example of the residuals after the PSF subtraction, 
for a resolved target, is presented in Fig. [3l upper panel. The resolved targets were then 
fitted with a two component model (PSF+galaxy). At high redshift, it becomes increas- 
ingly difficult to distinguish between bulge and exponential disk models from the luminosity 
distributions. We assumed that the host galaxies can be represented as elliptical galaxies 
with a de Vaucouleurs r^/^ profile. This choice is motivated by the strong evidence at low 
redshift for the predominance of bulge dominated hosts of quasars (e.g. Hamilton et al. 
2002; Dunlop et al. 2003; Pagani et al. 2003). If instead of a bulge-dominated model, we 
adopt an exponential disk law, all the objects can still be reasonably well fitted. The main 
difference is that the host galaxy luminosity become ~ 0.3 mag fainter for a disk model, but 
this does not significantly affect the results of this paper. 

In only two cases (2QZJ124029-0010 and 2QZJ222702-3205), the host galaxy detection 
was considered good enough to be able to firmly pin down the effective radius i?eff- For all 
the other objects in the sample, we op ted to set .Reff equa. 1 to 5 kpc in the modeling, following 



the indication of previous works (e.g.. lFalomo et al.ll2004l ). We note that, within our adopted 



cosmological framework, 5 kpc corresponds to 0.61 arcsec ai z = 2.5, comparable with the 
typical angular resolution of our data. In order to evaluate the effect on the host luminosity 
of assuming a more compact host galaxy, we re-performed the modelling of all the resolved 
targets assuming R^s = 2.5 kpc. With only two exceptions, the host galaxy magnitudes are 
consistent within 0.3 mag. No systematic offsets are found (the difference of average values 
= -0.03 mag). The two exceptions are Q2125-4432 and 2QZ23175 1-3147 (Am = 0.4 and 
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0.8 respectively), both of which have very faint host galaxies, close to the envelope of the 
unresolved targets, and therefore are likely to have larger uncertainties. Furthermore, the 
models of 2QZ231751-3147 have a highly elongated profile shape, suggesting that this object 
is peculiar. The number of unresolved targets remains unchanged irrespective of the assumed 
Res, since the resolved/unresolved classification only depends on the comparison between 
the observed light profile and the modelled PSF, independently of any other assumption on 
the shape of the host galaxy. 

Targets showing no residuals in the light profile after the PSF subtraction exceeding the 
PSF model uncertainty were considered unresolved. An example of an unresolved target is 
provided in Fig. [3l lower panel. In these cases, we computed the profiles expected assuming 
a zero-ellipticity host galaxy with a de Vaucouleurs profile and i?efr = 5 kpc, and various 
N/H ratios. Among them, we adopted as the upper limit on the host galaxy luminosity the 
model that best matched the PSF model uncertainty. 

5. Results 

In Fig. m we report the image of an example quasar, the best-fitting PSF, the residual 
after PSF subtraction, and the residual after the fit with both the PSF and the galaxy model. 
The fitted parameters are given in Table [2J We have been able to resolve 11 out of the 16 
quasars. In Table [2], we also report the nuclear and the host galaxy absolute magnitudes 
and the nuclear to host (N/H) luminosity ratio for each quasar. All objects with N/H<12 
are resolved, while no host galaxy is resolved in objects with N/H > 12, highlighting its key 
role in the detection of the host galaxy emission. 

To compare the properties of the quasar hosts at different redshifts, it is preferable to 
compare data probing the same rest-frame wavelengths. The Ks-band at redshifts 2 < z < 3 
corresponds to rest-frame ~ 5500 — 7000 A, closely matching the R-band. Therefore, in 
order to refer low and high redshift data to the same band (and to minimize color and 
fc-corrections), we transformed the observed magnitudes into absolute magnitudes in the 
R-band. To perform the color and /c-correction transformations, we assumed an elliptical 
galaxy template (Mannucci et al. 2001) for the host galaxy, and a composite quasar spectrum 
(Francis et al. 1991) for the nucleus. 

All the resolved quasars in our study have host galaxies with luminosity ranging between 
M(R) ~ —22.5 and M(R) ~ —25.5, i.e. corresponding to a range between ~ and 
~ M* — 3, where M=k(R)~ —22.5 is the R-band characteristic luminosity of the Schechter 
luminosity function for 2 < z < 2.5 inactive galaxies (Marchesini et al. 2007), which is ~ 1 
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mag brighter than the one observed in the local Universe (M*(R)~ —21.3; see, e.g., Gardner 
et al. 1997; Nakamura et al. 2003). The average R-band host luminosity of the resolved 
quasars is < M(R) host >= —24.2 ± 0.8. Following the approach of statistical survival 
analysis proposed by Feigelson & Nelson (1985), we have computed the average M(R)host 
including the upper limits of the host galaxy luminosities for the unresolved objects, and 
obtain < M(R) host >= —23.8 ± 0.5. For a more thorough discussion of this method, 
see Hyvonen et al. (2007a). For reference, the average nuclear luminosity of the quasars 
is < M(R) nucleus >= —25.5 ± 0.4, and the average nucleus-to-host luminosity ratio < 
N/H >= 6.7 ± 1.2. 

In Fig. [5] we compare the host galaxy luminosities derived in this paper with those avail- 
able from the literature for quasars in the same redshift range. All the published magnitudes 
are converted to rest-frame R-band assuming the elliptical template from Mannucci et al. 
(2001), and adopting the same cosmological framework. In the following, we compare our re- 
sults with each available sample in the literature. The targets in Ridgway et al. (2001) were 
selected for having very low quasar luminosity. Their two z ~ 2.7 quasars have M(R)host 
= —22.3 and —23.4 respectively, towards the faint end but consistent with the range of lu- 
minosities found in this paper. The lensed quasars in Peng et al. (2006) have host galaxies 
with luminosities very similar to those found in our sample, with < M(R) >= —24.2 ± 0.6 
(considering only the 2 < z < 3 targets). The same is true for the small number of resolved 
targets in the previous AO studies by our group (Falomo et al., 2005; 2008). Villforth, Heidt 
& Nilsson (2008) reported the marginal detection of the host galaxy of a single z = 2.75 
quasar. Its very low N/H ratio and spectral classification halfway between Type 1 and Type 
2 AGN suggest that this object is an obscured quasar. Due to this peculiarity, we excluded 
it from further analysis. Finally, Schramm et al. (2008) found very high host galaxy lumi- 
nosities, M(R)host = —25.8, —26.4 and —26.8, for three quasars at z = 2.643, z = 2.904 and 
z = 2.933, respectively. These luminosities are up to two magnitude brighter than in the 
other studies plotted in Fig. O This discrepancy may be connected with the fact that these 
quasars are among the most luminous at this redshift range. Motivated by this peculiarity, 
we have re-analyzed the Ks-band data of Schramm et al., retrieved from the ESO Archive. 
Our re-analysis shows that all these quasars are unresolved. In particular, no signal from the 
host galaxy was found at radii larger than ~ 2 arcsec (compare with Fig. 3 in Schramm et 
al., 2008). This is consistent with the presence of very bright nuclei in these quasars (about 
two magnitudes brighter than the brightest quasar in our sample). Furthermore, their obser- 
vations were performed in only modest seeing conditions (~0.7 arcsec FWHM). Therefore, 
we have excluded the Schramm et al. data from Fig. O and the following analysis. 

As noted previously, the nucleus-to-host luminosity ratio of the resolved targets in our 
sample ranges between 1 and 12. If quasars emit with a narrow range in Eddington ratio 
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and if the black hole mass is proportional to the luminosity of the host galaxy bulge (e.g., 
Magorrian et al., 1998), a correlation between the nuclear and the host galaxy magnitudes 
is to be expected. On the other hand, such a relation can fade due to several factors, both 
intrinsic (nuclear obscuration, beaming, and/or an intrinsic spread in the accretion rate) and 
extrinsic (e.g., lower level, non-quasar like nuclear activity may be relevant at low redshift; 
faint galaxies with very bright nuclei are difficult to resolve at high redshift), or can be 
an artifact of a narrowly sampled parameter space in the N/H ratio. Previous studies do 
not agree whether such a correlation exists (e.g., Dunlop et al., 2003; Pagani et al., 2003; 
Kotilainen et al., 2007). Fig. [6] shows the comparison between the rest-frame i?-band host 
and nuclear absolute magnitudes for all the available RQQ host galaxies in the 2 < z < 3 
range. No significant trend is observed between the host and the nuclear luminosity, neither 
considering only our sample nor comparing with all the available datasets, though we remark 
that the sampled parameter space in N/H strongly depends on the adopted observation 
technique, thus making a direct comparison between the various datasets difficult. We rather 
note that our new data, together with those published in Falomo et al. (2008) and most of 
the Ridgway et al. (2001) high redshift sample, fill a distinct region of the plot with respect 
to the data published in Peng et al. (2006). In other words, for a similar range in M(R)host, 
the quasars in Peng et al. have ~ 3 mag fainter nuclei, and their N/H< 1, suggesting that 
the sample of Peng et al. is dominated by lower luminosity objects with non-quasar level, 
e.g. Seyfert type, nuclear activity. 



6. Discussion 

In Fig. [7] we extend the plot shown in Fig. \5\ to include all published RQQ host 
galaxy magnitudes at 0.5 < z < 2 (Kukula et al. 2001; Ridgway et al. 2001; Falomo et 
al. 2004; Kotilainen et al. 2007; Hyvonen et al. 2007a), and HST observations at z < 1 
(Bahcall et al., 1997; Hooper et al., 1997; Boyce et al., 1998; Kukula et al., 2001; Hamilton 
et al, 2002; Dunlop et al, 2003; Floyd et al., 2004; Labita et al, 2006). In order to treat 
these literature data homogeneously, we transformed the published apparent magnitudes 
to M(R)host following the procedure described above (/c-correction, cosmology and color 
correction). Note that Groom et al. (2004) resolved the host galaxy of one quasar out of 
nine observed. Since they do not report upper limits on the magnitudes of the unresolved 
quasars, it is not possible to assess whether their unresolved hosts are consistent with the 
trend shown in Fig. [71 Because of this, the contribution from the single detected host galaxy 
by Groom et al. is negligible in this context and was omitted. 

We are now in a position to compare the observed general trend of the luminosity 
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of quasar host galaxies as a function of cosmic epoch with the expectations of theoretical 
models of galaxy formation. Fig. [7] shows that, despite the presence of a considerable scatter, 
a general trend is apparent, with the host galaxy luminosity of RQQs increasing by ~ 1.5 mag 
from the present epoch up to z ~ 3. We stress that this trend is statistically significant. For 
example, considering data binned as median averages, in order to be unaffected by the true 
(unknown) magnitude of the unresolved quasar hosts (see Figure [8]), the linear best-fit of the 
relation is: M(R)host= -0.56±0.1 x z - 22.5±0.3, with = 1.3 (where the uncertainties are 
3a errors from the maps and the is normalized to the number of degrees of freedom ). 
On the other hand, a fit with constant luminosity implies M(R)host = -23.2, with = 31.4. 
Therefore, the trend in Fig. [7] can not be reproduced with the host luminosity remaining 
constant with redshift. 

Most semi-analytical hierarchical models predict very few old, massive galaxies at high 
redshift (z > 1) because in these models, large structures (massive galaxies) preferentially 
form at late epoch by continuous merging of smaller galaxies (e.g. Cole et al. 2000; Croton 
et al. 2006). This trend clearly disagrees with the observations presented here (Fig. [7j), and 
with the discovery of a substantial population of evolved, massive galaxies at z > 1.5 (e.g. 
Daddi et al. 2005). Recent hierarchical models include feedback effects from supernovae 
and AGN to disentangle the baryon evolution from the hierarchical assembly of dark matter 
structures (e.g., Di Matteo et al., 2005; Bower et al. 2006; De Lucia et al., 2006; see Elhs, 
2008 for a review on this topic ). AGN feedback acts expelling gas from the galaxies and 
thus quenching star formation (ITremonti et al.l 120071 : iBundy et al.ll2008l ). This effect is of 



fundamental importance in massive galaxies, where AGN activity is preferentially observed 
(Kauffmann et al., 2003; Decarh et al., 2007; Gallo et al., 2008). 

In order to probe whether inactive galaxies and quasar host galaxies have different star 
formation histories, we overplot in Fig. [H]the trend observed for the characteristic luminosity 

MUDS 

function of redshift, as derived from the UKIDSS Ultra Deep Survey (Lawrence 
et al., 2007). The evolution of the luminosity function (Cirasuolo et al. 2007, 2008) is of the 
form: 

/ N (0.47±0.2) 

logMUDSAB(^) = -22.26 - (1) 

^ * ^ ' Vl-78±0.15y ^ ' 

We transformed the Ks"^^ magnitudes to the Ks-band assuming the correction used in the 
GOODS-ISAAC survey (Grazian et al. 2006): Ks^^=Ks + 1.895. Then we applied filter and 
^-correction as described above to compute the rest-frame R-band absolute magnitude. The 
result is shown in Fig. [HI On average, the available data on RQQ host galaxy luminosities 
lie between M^^s ^j^fi M^^^ — 1 at all redshifts. This reinforces the leading assumption 
that we make that quasars are harboured in luminous, massive galaxies. Furthermore, the 
dependence of M^, on z for quiescent galaxies (Cirasuolo et al. 2007, 2008) remarkably closely 
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matches the trend observed in the M(R)host - z distribution for quasar host galaxies, and 
is consistent with the findings of other studies on the luminosity function of galaxies at high 
redshift (e.g., Marchesini et al. 2007). This dependence supports the idea that quasar hosts 
and massive inactive galaxies share a similar star formation history. 

There is increasing evidence for a mass dependence of the star formation history of 
galaxies (e.g., Gavazzi et al., 1996; 2002; Cimatti et al., 2004; Treu et al., 2005a; 2005b; 
Thomas et al., 2005), in the sense that the more massive elliptical galaxies formed their stars 
in relatively shorter bursts of intense star formation and at higher redshift (z ~ 2.5 — 5), 
compared to less massive galaxies. Recent merger events, if and when they occur, do not 
significantly affect the content of the stellar population of massive galaxies. On the other 
hand, the formation of lower luminosity galaxies is shifted toward lower redshifts, having 
significant star formation present at all epochs (see e.g. Renzini 2006, Scarlata et al. 2007, 
Vergani et al. 2008). 

We stress that while the case of the evolution of a single burst population formed at high 
redshift provides an adequate explanation for our results, we can not rule out other models, 
accounting for the accretion of the galaxy mass due to mergers and different stellar history 
recipes. While a naive model with galaxies experiencing substantial mass accretion is not 
consistent with the data, a more complex model in which episodic star formation counter- 
weights the mass accretion of galaxies is formally acceptable. Nevertheless, we prefer an 
interpretation in which such a coincidence is not necessary. A more detailed modelling of 
the evolution and the star formation history of the host galaxies, that takes into account 
the role of the active nucleus is required. However, that is beyond the scope of this paper. 
Finally, we note that the assumption of a different SED of the stellar population does not 
affect significantly our results, since the rest-frame R-Ks color addressed in this paper is 
fairly insensitive to the age of the stellar population. 



7. Summciry and conclusions 

We have presented homogeneous high resolution NIR images for a sample of 16 low 
luminosity RQQs in the redshift range 2 < z < 3, to characterize the properties and the 
cosmological evolution of their host galaxies, in conjunction with data at lower redshift. The 
host galaxy was resolved in eleven quasars out of the 16, while the remaining were unresolved. 
The RQQs in our sample have < M(R) host >= -23.8±0.5, < M(R) nucleus >= -25.5±0.4 
and<N/H >= 6.7 ±1.2. 

Comparing our new data with literature data at 2 < 2; < 3 and at lower redshift. 
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we found that the quasar host galaxies follow the trend in luminosity of massive inactive 
ellipticals (with M(R) ranging between M^, and — 3 at 2 < ^ < 3, or — 1 and — 4 at 
z=0) undergoing passive evolution. This represents a fundamental constraint for models of 
galaxy formation and evolution and throws new light on our understanding of the interplay 
between the black holes, nuclear activity and the host galaxy evolution. 

To definitely pin down the host galaxy assembly history, a large observational effort 
in increasing the sample of resolved host galaxies at 2; > 2 — 3 is required, in order to 
increase the statistical significance of the present results and to sample a wider parameter 
space, whether in the N/H ratio, in the host galaxy luminosity, or in the level of quasar 
radio emission. Further observations with very high S/N and a very narrow, extremely 
reliable PSF, eventually thanks to laser guide star assisted AO in NIR will play a crucial 
role. Improved knowledge about the colours and environments of quasar host galaxies will 
also be needed to test the predictions of galaxy formation models. Finally, in a forthcoming 
paper (Decarli et al., in preparation), we will compare the luminosity of the host galaxies of 
high redshift quasars with their black hole masses, as derived from the spectroscopy of their 
broad emission lines. This will enable a direct test of the evolution of the black hole mass - 
host galaxy relations up to 2; = 3. 
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ison with inactive galaxies, and for the anonymous referee for constructive criticism. This 
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Fig. 1. — The rest-frame V-band absolute magnitude, as derived from the observed apparent 
V-band magnitude and assuming the Francis et al. (1991) template for /c-correction, as a 
function of redshift for our sample objects (filled circles) and for similar objects in the 
literature (see the inserted symbol caption for details). Solid lines show the cut in absolute 
magnitude that we adopted in the target selection. The contours (dotted lines) represent 
the distribution of all the quasars in the quasar catalogue of Veron-Cetty & Veron (2006), to 
which we refer to for details. The contour levels are logarithmically spaced, and reflect the 
number of quasars within each contour. Note that the absolute magnitudes of the targets 
from Peng et al. (2006) are corrected for the magnification due to gravitational lensing. 
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Fig. 2. — An example of the PSF modeling, for the frame of 2QZJ124029-0010. The average 
and standard deviation of the FWHM and ellipticity of the stars are listed, a) the contour 
plot of a star in a 40 x 40 pixel (6x6 arcsec) box. b) the same star in a 3D plot, c) the 
positions of the stars (filled circles) used in the PSF modeling within the frame. The position 
of the target is also shown as an empty diamond, d) the fitted PSF (solid line), compared 
to the observed star profiles (filled circles), e) the residuals after the PSF modeling. The 
shaded area represents the estimated uncertainty in the PSF model, derived as described in 
the text. 
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Fig. 3. — Tlie residuals after the PSF fitting of tlie target (blue, empty diamonds) compared 
to the uncertainties in the PSF model (shaded area, see also panel e in Fig. [2]). In the 
upper and lower panel, we show an example of a resolved target (2QZJ124029-0030) and an 
unresolved target (SDSSJ16187-0043), respectively. The green, dot-dashed lines refer to the 
simulated residual profiles for a quasar with a host galaxy following a de Vaucouleurs profile 
with i?efr = 5 kpc and various N/H ratios. 
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Fig. 4. — The central 4x4 arcsec region surrounding the quasar 2QZJ124029-0010. In 
the left and middle panels, from top to bottom: (a) the original image, (b) the PSF model, 
(c) the host galaxy (PSF model subtracted from the observed profile), and (d) residuals of 
the fit. On the right, the top panel shows the observed radial profiles of the quasars (open 
diamonds), superimposed to the PSF model (blue dotted line) and an elliptical galaxy model 
convolved with its PSF (green dashed line). The red solid line shows the composite fit. The 
second panel shows the radial profiles of the stars compared to the PSF model (red solid 
line). The third panel shows the radial profile of the host galaxy after the PSF subtraction, 
while the bottom two panels show the residuals and the distribution of the fit. 
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Fig. 5. — The host galaxy absolute magnitudes as a function of redshift for the targets in our 
work (filled circles) compared with all the available literature data in the 2 < 2; < 3 redshift 
range. The arrows represent the upper limits of the host luminosity for the unresolved 
quasars in our study. The dash-dotted and dotted lines show, for comparison, the values of 
M^, and M^, — 3, where the characteristic R-band magnitude of the luminosity function 
of 2 < z < 2.5 inactive galaxies (Marchesini et al. 2007). 
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Fig. 6. — The absolute magnitude of the RQQ host galaxies compared to that of the nucleus 
in our sample (filled circles) and in other samples from the literature. The diagonal lines 
represent the loci of constant ratio between host and nuclear emission at N/H=l (solid line) 
and N/H=10 (dashed line). For a similar range in M(R)host, tha data seem to fill two 
different regions of the plot, one with —24 > M(R)nucleus > —26 (data from Ridgway et 
al., 2001, Falomo et al., 2008 and this work) and one with -20 > M(R)nucleus > -23 
(dominated by the sample of Peng et al., 2006). Note that having N/H<1, objects in the 
sample of Peng et al. should strictly be classified as lower luminosity AGNs (e.g. Seyferts) 
instead of quasars. 
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Fig. 7. — The redshift distribution of the host galaxy luminosity of radio-quiet quasars 
in this paper (filled circles) compared with all the available literature data. The arrows 
represent the upper limits of the host luminosity for the unresolved quasars. An overall 
~ 1.5 mag increase of the average M(R)host is observed from = to 2 = 3. 
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Fig. 8. — The redshift distribution of the host galaxy luminosity of radio-quiet quasars 
plotted in Fig. [TJ here binned as median averages at z = 0.5, 1.5 and 2.5. Uncertainties are 
the quartile values of the data distribution in each bin. The upper panel shows a histogram 
of the number of targets per bin. The shaded area represents the evolution of the knee of 
inactive galaxy luminosity function from UKIDSS Ultra Deep Survey (Cirasuolo et al. 2007, 
2008), /c-corrected assuming the SED of a single-burst stellar population with ^^burst = 5. 
The area extends from M^^^ to M^°^ — 1. We note that the R-Ks color is poorly sensitive 
to the age of the stellar population: The difference between the adopted correction and a 
fixed ~ 2.65 color is < 0.15 mag for z ranging between and 3. 
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Table 1. Journal of the observations 



Quasar 




z 


ya 


Date (2006) 




Seeing 








ZP'i 




sky 


mag'^ 








[mag] 






[arcsec] 






[mag] 




[mag/arcscc^] 




2QZJ124029-0010 


2. 


,030 


19.76 


11- Jul 





.52 


±0, 


.02 


2 


23. 


,69±0, 


.04 


12. 


,85 


2QZJ133136-0002 


2. 


,710 


20.74 


10- Jul 





.52 


±0, 


,01 


4 


24, 


,09±0, 


,06 


13. 


,30 


2QZJ143220-0215 


2. 


,476 


20.38 


17-May 





.42 


±0, 


,01 


5 


24, 


,14±0, 


,08 


13, 


,39 


2QZJ144022-0122 


2. 


,244 


20.02 


21-May, 10- Jul 





.46 


±0, 


,05 


1 


24, 


,02±0, 


,01 


13, 


,53 


SDSSJ16187-0043 


2. 


,068 


19.77 


13- Apr 





.45 


±0, 


,00 


4 


24, 


,00±0, 


,08 


13, 


,48 


Q2125-4432 


2. 


,503 


20.39 


13-May 





.56 


±0, 


,01 


3 


24. 


,07±0, 


,02 


13. 


,00 


Q 2126-1148A 


2. 


,188 


20.00 


13-May 





.45 


±0, 


.01 


5 


24. 


,12±0, 


.02 


13. 


,26 


2QZJ215539-3026 


2. 


,593 


20.44 


17-May 





.43 


±0, 


.02 


6 


24. 


,10±0, 


.08 


13. 


,04 


2QZJ221139-3132 


2. 


,391 


20.40 


17-May 





.47 


±0, 


.02 


5 


24. 


,14±0, 


.09 


13. 


,22 


Q2225-403 


2. 


,410 


20.20 


18-May, 11-Jun 





.46 


±0, 


.05 


1 


24. 


,05±0, 


.04 


13. 


,34 


2QZJ222702-3205 


2. 


,177 


20.13 


12-Jun 





.45 


±0, 


,00 


2 


24, 


,09±0, 


,00 


13, 


,42 


2QZJ223048-2954 


2, 


,652 


20.56 


12-Jun 





.51 


±0, 


,02 


4 


24, 


,14±0, 


.14 


13. 


,44 


2QZJ225950-3206 


2. 


,225 


19.72 


17- Jul 





.56 


±0, 


,01 


2 


23, 


,95±0, 


,10 


13, 


,24 


2QZJ231751-3147 


2. 


,628 


20.58 


18- Jul, 2-Aug 





.56 


±0, 


.01 


6 


24. 


,10±0, 


,03 


13. 


,22 


2QZJ232755-3154 


2. 


,737 


20.73 


14-Jul, 17- Jul 





.40 


±0, 


.01 


4 


23. 


,99±0, 


.04 


13. 


,11 


2QZJ233451-2929 


2. 


,669 


20.72 


19-Jun, 3-Aug 





.46 


±0, 


.02 


2 


24. 


,16±0, 


.12 


13. 


,46 



^Quasar V-band apparent magnitudes from Veron-Cetty & Veron (2006). 

^The average and rms FWHM in arcsec of all stars in the frame. When only 1 non-saturated star 
was present, from the comparison of the PSF with resolved targets we conservatively estimated a 
~ 0.05 arcsec uncertainty in the seeing. 

■^Number of stars used in the seeing estimate. 

'^Zero Point from 2MASS magnitude of field stars. 

^Sky surface brightness. 
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Table 2. Properties of the quasar host galaxies. 



Quasar 



m(Ks)^ m(Ks)^ M(R)^ M{Ry N/ff 
nuc host nuc host 



2QZJ124029-0010 


2 


.030 


17, 


.2 


17.4 


-25. 


.9 




■25.6 


1.3 


2QZJ133136-0002 


2 


.710 


18, 


.4 


> 21.3 


-25. 


.6 


> 


-22.6 


>16.0 




9 




1 s 

lo. 




.o 


9n ^ 


9^ 
-zo. 


9 

,z 




9*? 1 


u. / 


2QZJ144022-0122 


2 


.244 


17, 


.7 


18.7 


-25. 


,5 




-24.6 


2.2 


SDSSJ16187-0043 


2 


.068 


17, 


.0 


> 20.4 


-26. 


,0 


> 


-22.6 


>23.1 


Q2125-4432 


2, 


.503 


18, 


.2 


20.7 


-25. 


,4 




-22.9 


9.5 


Q 2126-1148A 


2, 


.188 


17, 


.2 


20.1 


-25. 


,8 




-23.1 


12.0 


2QZJ215539-3026 


2, 


.593 


18, 


.3 


19.5 


-25. 


.6 




■24.2 


3.4 


2QZJ221139-3132 


2 


.391 


17, 


.6 


> 20.5 


-25. 


.7 


> 


-23.0 


>12.2 


Q2225-403 


2 


.410 


17, 


.4 


18.5 


-25. 


.9 




■25.0 


2.4 


2QZJ222702-3205 


2 


.177 


18, 


.2 


18.5 


-24. 


.8 




■24.7 


1.2 


2QZJ223048-2954 


2 


.652 


18, 


.6 


> 21.5 


-25. 


,4 


> 


-22.4 


>15.8 


2QZJ225950-3206 


2 


.225 


17, 


.5 


18.6 


-25. 


,6 




-24.7 


2.4 


2QZJ231751-3147 


2 


.628 


19, 


.2 


20.4 


-24. 


,8 




-23.4 


3.5 


2QZJ232755-3154 


2 


.737 


18, 


.6 


19.5 


-25. 


,5 




-24.5 


2.5 


2QZJ233451-2929 


2, 


.669 


18, 


.5 


> 21.9 


-25. 


,5 


> 


-22.0 


>27.3 



'^Apparent magnitudes of the nuclei in the observed Ks-band. 

''Apparent magnitudes of the host galaxies in the observed Ks-band. 

'^Absolute magnitudes of the nuclei in the R-band, /c-corrected as- 
suming the Francis et al. (1991) template; no correction for galactic 
extinction is applied. 

'^Absolute magnitudes of the host galaxy in the R-band, /c-corrected 
assuming the eUiptical galaxy template by Mannucci et al. (2001); no 
correction for galactic extinction is applied. 

^The N/H ratio, referred to the absolute R magnitudes. 



